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Abstract The general features of the cosmic-ray spectrum have been known for a 
long time. Although the basic approaches to understanding cosmic-ray 
propagation and acceleration have also been well understood for many 
years, there are several questions of great interest that motivate the cur- 
rent intense experimental activity in the field. If the energy-dependence 
of the secondary to primary ratio of galactic cosmic rays is as steep as 
observed, why is the flux of PeV particles so nearly isotropic? Can all 
antiprotons and positrons be explained as secondaries or is there some 
contribution from exotic sources? What is the maximum energy of cos- 
mic accelerators? Is the "knee" of the cosmic-ray spectrum an effect 
of propagation or does it perhaps reflect the upper limit of galactic ac- 
celeration processes? Are gamma-ray burst sources (GRBs) and/or ac- 
tive galactic nuclei (AGN) accelerators of ultra-high-energy cosmic rays 
(UHECR) as well as sources of high-energy photons? Are GRBs and/or 
AGNs also sources of high-energy neutrinos? If there are indeed par- 
ticles with energies greater than the cutoff expected from propagation 
through the microwave background radiation, what are their sources? 
The purpose of this lecture is to introduce the main topics of the School 
and to relate the theoretical questions to the experiments that can an- 
swer them. 

Keywords: High Energy Cosmic Rays 

1. Introduction 

It is appropriate to begin this lecture with a diagram from the review 
of Shapiro & Silberberg, 1970, which compares the abundances of ele- 
ments in the cosmic radiation with solar system abundances. This classic 
measurement is one of the foundations of cosmic-ray physics. The ele- 
ments lithium, beryllium and boron are quite abundant among cosmic 
rays even though they constitute only a tiny fraction of the material in 
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Figure 1. Comparison of cosmic-ray abundances with abundances of nuclei in the 
solar system (Shapiro & Silberberg, 1970). 



the solar system and the interstellar medium. This fact is understood 
largely as the result of spallation of the "primary" nuclei carbon and 
oxygen during their propagation in the interstellar medium. The idea is 
that the protons and nuclei are accelerated from the interstellar medium 
and/or from the gas in or near their sources. The composition of the 
"primary" nuclei (defined as those initially accelerated and injected into 
the interstellar medium with high energy) reflects the combination of 
nuclei present in the material that is accelerated, which contains negligi- 
ble amounts of "secondary" elements such as Li, Be, B. These secondary 
nuclei are fragmentation products of heavier primary nuclei. 

In this picture, the amount of secondary nuclei is a measure of the 
characteristic time for propagation of cosmic rays before they escape 
from the galaxy into inter-galactic space. A simplified version of the 
diffusion equation that relates the observed abundances and spectra to 
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initial values is 
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. ... Q l (E)-(pcn H a i + —) N z (E)+pcn H Y, (J k^iN k (E). (1) 

Tesc(E) \ mJ 

Here Ni{E) is the spatial density of cosmic-ray nuclei of mass i, and 
iih is the number density of target nuclei (mostly hydrogen) in the in- 
terstellar medium, Qi(E) is the number of primary nuclei of type i ac- 
celerated per cm 3 per second, and <7j and au^i are respectively the total 
and partial cross sections for interactions of cosmic-ray nuclei with the 
gas in the interstellar medium. The second term on the r.h.s. of Eq. 1 
represents losses due to interactions with cross section <jj and decay for 
unstable nuclei with lifetime r«. The energy per nucleon, E, remains 
constant to a good approximation in the transition from parent nuclei 
to nuclear spallation products, which move with velocity (3c and Lorentz 
factor 7 = E/m p . 

A crude first estimate of the characteristic diffusion time r esc can 
be made by neglecting propagation losses for a primary nucleus P and 
assuming that Qs = for a secondary nucleus S. If we also neglect 
collision losses by the secondary nucleus after it is produced, then the 
solution of Eq. 1 is 

n H T esc = -r— ^ -p. (2) 

Since the density in the disk of the galaxy is of order one particle per 
cm 3 and the typical partial cross section for a light nucleus will be of 
order 100 mb, the scale for the characteristic time is ~ 10 7 years. The 
full analysis requires self-consistent solution of the coupled equations for 
all species accounting for all loss terms, as described in John Wefel's 
lecture in this volume (Wefel, 2005). A nice overview of propagation 
models is given by Jones et al, 2001. 

The thickness of the disk of the galaxy is of order 300 pc = 1000 light 
years, which is much shorter than the characteristic propagation time 
of 10 million years. The explanation is that the charged particles are 
trapped in the turbulent magnetized plasma of the interstellar medium 
and only diffuse slowly away from the disk, which is assumed to be 
where the sources are located. Measurements of the ratio of unstable to 
stable secondary nuclei (especially 10 Be/ 9 Be) are used to determine T esc 
independently of the product nn T esc and hence to constrain further the 
models of cosmic-ray propagation. 

Another important fact is that the ratio N$/Np is observed to de- 
crease with energy. From Eq. 2 this implies that r esc also decreases. 
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Simple power-law fits to ratios like B/C give 

7~esc 0^ E i 



(3) 



with 5 ~ 0.6. This behavior has important consequences for the source 
spectrum. To see this, consider an abundant, light primary nucleus such 
as hydrogen or helium. They are sufficiently abundant so that feed- 
down from heavier nuclei can be neglected and their cross sections are 
small enough so that energy losses in the interstellar medium can also 
be neglected in a first approximation. Then Eq. 1 reduces to 



The local energy-density spectrum of cosmic rays is related to the ob- 
served flux (^(particles per cm 2 per GeV per second per steradian) by 



Since 4>(E) ps E~ a with a ps 2.7 the inference is that the cosmic accel- 
erators are characterized by a power law with spectral index as ~ 2.1. 
This value is close to the spectral index for first order acceleration by 
strong shocks in the test-particle approximation (Ostrowski, 2005). 

If the time for diffusion out of the galaxy continues to be described 
by 5 ps 0.6 to very high energy there is a problem: as cr esc decreases 
and approaches galactic scales, the cosmic-ray fluxes should become sig- 
nificantly anisotropic, which is not observed. One possibility is that 
the observed energy-dependence at low energy is due to a combination 
of "reacceleration" (Seo & Ptuskin, 1994, Heinbach & Simon, 1995) by 
weak shocks in the interstellar medium after initial acceleration by strong 
shocks. In such models, the high-energy behavior of diffusion is typically 
described by a slower energy dependence with 5 ps 0.3. If so, the source 
spectral index would be steeper, approximately a s = 2.4. In any case, in 
more realistic non-linear treatments of acceleration by strong shocks, the 
spectrum has some curvature, being steeper at low energy and harder at 
the high energy end of the spectrum (Berezhko & Ellison, 1999). What 
we observe may be some kind of average over many sources, each of 
which is somewhat different in shape and maximum energy. 

The assumption underlying the discussion above is that the sources 
accelerating cosmic rays are in the disk of the galaxy and that the energy 
density in cosmic rays observed locally is typical of other regions of 
the galactic disk. If so, the total power Pcr required to maintain the 
cosmic radiation in equilibrium may be obtained by integrating Eq. 4 
over energy and space. The result is 



Pcr = / d 3 x / Q(E)dE = V G — <j>{E) / r esc {E) dE. (6) 



Q(E) = N(E)/ Tesc ~ N(E) x E s . 



(4) 



N(E) = —<j>(E). 



(5) 
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Using the observed spectrum and the value of r esc explained above, one 
finds numerically 



The kinetic energy of an expanding supernova remnant is initially of 
order 10 51 erg, excluding the neutrinos, which carry away most of the 
energy released but do not disturb the interstellar medium. There are 
perhaps 3 supernovae per century, which gives Psn ~ 10 42 erg/s as an 
estimate of the power which is dissipated in the interstellar medium by 
means of strong shocks driven by supernova ejecta. These are just the 
ingredients needed for acceleration of galactic cosmic rays. I return to 
the subject of cosmic-ray acceleration in §3 below. 

2. Secondary cosmic-rays 

Depending on the context, the term secondary cosmic rays can refer 
either to particles produced by interactions of primary cosmic rays with 
the interstellar gas or to particles produced by interactions of cosmic rays 
in the Earth's atmosphere. The production mechanisms are similar, and 
there are some common features. 

We have already seen one example, the secondary nuclei produced by 
occasional interactions of primary cosmic-ray nuclei during their propa- 
gation in the interstellar medium. In that case, to a good approximation, 
the energy per nucleon of the secondary nucleus is the same as that of 
the parent nucleus. The reason is that the nuclear fragments are only 
spectators to any production of secondary pions that may occur in the 
collisions. In general, however, the production term (last term of Eq. 1) 
will involve an integration over the energy of the parent particle. The 
main process to consider is production of pions by interaction of protons 
with a target nucleus. The production spectrum of the pions is 



Here X p = Am p /o™ x is the interaction length of protons in a medium 
consisting of nuclei of mass A, and dX = pdi is the differential element 
of mass traversed in distance d£ in a medium of density p. It is often a 
useful approximation at high energy (i.e. energy 3> particle masses) to 
assume a scaling form for the dimensionless production spectrum: 



P CR ~ 10 41 erg/s. 



(7) 






dn^ (E n , E p ) 
dE w 



(9) 
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The scaling variable is £ = E n /E p , and Eq. 8 becomes 

dx f\ ( E A <M0 d ^ dX ^ IT <> 7 



(10) 

The last step on the r.h.s. of Eq. 10 follows when the parent spectrum is 
a power law in energy {(j) p {E) = K E~ a ). In that case, in the high-energy 
scaling approximation 

d(p n (E w ) -» ^-Zp-n x (f> p (E n ), (11) 

i.e. the energy spectrum of secondaries has the same power behavior as 
the primaries scaled down by a factor 

W) = J^V-^df. (12) 

The spectrum- weighted moment Z p ^{a) depends both on the physics 
of production of the secondary pion and on the value of the differential 
spectral index a. 



2.1 Galactic secondaries 

2.1.1 Diffuse gamma-rays. Gamma-ray emission from the 
disk of the galaxy is a powerful probe of the model of cosmic-ray origin 
and propagation as well as of the structure of the galaxy. Unlike charged 
cosmic rays, secondary photons propagate in straight lines. Since the 
galaxy is transparent for 7-rays of most energies, it is possible to search 
for concentrations of primary cosmic-ray activity from the map of the 
7-ray sky after subtracting point sources. For example, if cosmic-ray 
acceleration is correlated with regions of higher density such as star- 
forming regions where supernovae are more frequent, then one would 
expect a quadratic enhancement of secondary production because of the 
spatial correlation between primary flux and target density. 

The baseline calculation is to assume that the intensity observed lo- 
cally at Earth is representative of the distribution everywhere in the disk 
of the galaxy. One can then look for interesting variation superimposed 
on this baseline. Following the analysis of Eqs. [9,10,11,12], the aver- 
age number of neutral pions produced per GeV per unit volume in the 
interstellar medium is 



q n = 4irn H a p n £ Z^o (f> p (E n ). (13) 
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Next this expression has to be convolved with the distribution of pho- 
tons produced in tt° — > 7 7. In the rest frame of the parent pion, each 
photon has E 1 = = 70 MeV. The angular distribution is isotropic, 

so 

dn 7 1 dn 7 1 

dD* ~ 2^ d cos 6* 2tt' ^ 
where 6* is the polar angle of the photon along the direction of motion 
of the parent pion but evaluated in the rest frame of the pion. For decay 
in flight of a pion with Lorentz factor 7 and velocity /3c, the energy of 
each of the resulting photons is 

£ 7 = 7^(1 + 0cos9*) (15) 

with cos 0* = — cos 62 for the two photons. Changing variables in Eq. 14 
then gives 

dE^ P^rrin 

For E n > E^ ^> the convolution of the distribution (16) with the 

production spectrum of neutral pions (13) gives 

g 7 (£ 7 ) « 4vr n H a™ 1 x ^ Z p ^o(a) <j> p {E^). (17) 

Numerically, in the approximation of uniform cosmic-ray density and 
uniform gas density in the interstellar medium, the observed gamma-ray 
flux would be 

( E "t) ~ o m-6 ( ' n H\ ( rmax(b,(f>e) 



^gi „ 3 x 1(T 6 ^ 'm°*V»>W . (18) 
(p p {E~f) \cm 6 ) \ lkpc J 

Here r max {b,4>i) is the distance in a direction {b,<f>i} to the effective 
edge of the galactic disk, where b and <pe are galactic latitude and lon- 
gitude. The effective distance is defined with respect to an equivalent 
disk of uniform density. Eq. 18 compares well in order of magnitude 
with the measured intensity of GeV photons by Egret (Hunter et a/., 
1997). The derivation of Eq. 18 is grossly oversimplified compared to the 
actual model calculation made in the paper of Hunter et ah, 1997 of the 
diffuse, galactic gamma-radiation. The reader is urged to consult that 
paper to understand the impressive level of detail at which the data are 
understood. It is also interesting to compare Eq. 18 to the TeV diffuse 
flux measured by Milagro (Goodman, 2005). 

The implication of Eq. 17 is that for E^ 3> 70 MeV the diffuse gamma- 
ray spectrum should have the same power law behavior as the proton 
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spectrum, a « 2.7. What is observed, however, is that the spectrum of 
gamma-rays from the inner galaxy is harder than this, having a power- 
law behavior of approximately E~ 2A (Hunter et al, 1997). This is 
currently not fully understood. One possibility is that the cosmic-ray 
spectrum producing the gamma rays is harder than observed locally near 
Earth (Hunter et al, 1997). 

Cosmic-ray electrons also contribute to the diffuse gamma-radiation 
by bremsstrahlung and by inverse Compton scattering. Fitting the ob- 
served spectrum requires a complete model of propagation that includes 
all contributions (Hunter et al, 1997). The distinguishing feature of 
7r°-decay photons is a kinematic peak at E y = m n /2. The origin of this 
feature can be seen in Eq. 15 from which the limits on E y for any given 
Lorentz factor of the parent pion are 



The distribution dn^/dEj is flat between these limits for each 7 and 
is always centered around ln(m 7r /2) when plotted as a function ln(.E 7 ). 
The individual contributions for parent pions of various energies always 
overlap at ln(.E 7 ) = m n /2, so the full distribution from any spectrum 
always peaks at this value (Stecker, 1971). 

2.1.2 Antiprotons and positrons. Antiprotons and positrons 
are of special interest because an excess over what is expected from pro- 
duction by protons during propagation could reflect an exotic process 
such as evaporation of primordial black holes or decay of exotic relic 
particles (Bottino et al., 1998). At a more practical level, they are im- 
portant because they are secondaries of the dominant proton component 
of the cosmic radiation. As a consequence their spectra and abundances 
provide an independent constraint on models of cosmic-ray propagation 
( Moskalenko et al, 1998). 

Secondary antiprotons have a kinematic feature analogous to that in 
7r°-decay gamma rays but at a higher energy related to the nucleon mass. 
In this case the feature is related to the high threshold for production of 
a nucleon-antinucleon pair in a proton-proton collision. This kinematic 
feature is observed in the data (Orito et al., 2000), and suggests that 
an exotic component of antiprotons is not required. Antiproton fluxes 
are consistent with the basic model of cosmic-ray propagation described 
in the Introduction. 

Positrons are produced in the chain 




P — > 7T 
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Secondary electrons are produced in the charge conjugate process, but 
their number in the GeV range is an order of magnitude lower than 
primary electrons (i.e. electrons accelerated as cosmic rays). Because 
of radiative processes, the spectra of positrons and electrons are more 
complex to interpret than high-energy secondary 7-rays (Moskalenko 
k, Strong, 1998). The measured intensity of positrons appears to be 
consistent with secondary origin (DuVernois et al, 2001). 

2.1.3 7-rays and u from young supernova remnants. The 

same equations that govern production of secondaries in the interstellar 
medium also apply to production in gas concentrations near the sources. 
For example, a supernova exploding into a dense region of the interstellar 
medium (Berezhko & Volk, 2000) or into the gas generated by the strong 
pre-supernova wind of a massive progenitor star (Berezhko, Piihlhofer 
& Volk, 2003) would produce secondary photons that could show up as 
point sources. 

Indeed, for many years, observation of 7T°-decay 7-rays from the vicin- 
ity of shocks around young supernova remnants (SNR) has been consid- 
ered a crucial test of the supernova model of cosmic-ray origin (Drury 
et al., 1994; Buckley et al., 1998). Note, however, that a sufficiently 
dense target is required. Moreover, it is difficult to distinguish photons 
from 7T° decay from photons originating in radiative processes of elec- 
trons ( Gaisser, Protheroe & Stanev, 1998). There are two signatures: 
at low energy, observation of a shoulder reflecting the 7r° peak at 70 
MeV would be conclusive. At higher energy one has to depend on the 
hardness and shape of the spectrum for evidence of hadronic origin of 
the photons. (See Drury et al, 2001 for a current review.) 

Observation of high-energy neutrinos would be strong evidence for 
acceleration of a primary beam of nucleons because such neutrinos are 
produced in hadronic interactions. Expected fluxes are low (Gaisser, 
Halzen Sz Stanev, 1995), so large detectors are needed (Montaruli, 2003, 
Migneco, 2005). 

2.2 Atmospheric secondaries 

Production of secondary cosmic rays and 7-rays in the interstellar 
medium generally involves less than one interaction per primary. In the 
language of accelerators, this is the thin-target regime. In contrast, the 
depth of the atmosphere is more than ten hadronic interaction lengths, 
so we have a thick target to deal with. The relevant cascade equation is 
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dNi(E,X) (Nj(E,X) | Nj(E,Xy . (J()) 



where 



ax V Ai(s) di{E) 

^ [°° N k (E',X)F k ^ i (E,E') AW/ 

+ \h -x^o 1 



1 dou^i 
?k^i - — E . 
Ok dE 



The equation describes the longitudinal development of the components 
of the atmospheric cascade in terms of slant depth (dX = p d£) along 
the direction of the cascade. 

The loss terms on the r.h.s. of Eq. 19 represent interactions and decay, 
in analogy to Eq. 1. Here 

di = pjCTi = p^-^r (20) 

is the Lorentz dilated decay length of particle i in g/cm 3 . The expression 
Aj = di defines a critical energy below which decay is more important 
than re-interaction. Because the density of the atmosphere varies with 
altitude, it is conventional to define the critical energy at the depth of 
cascade maximum (Gaisser, 1990). For pions the critical energy in the 
terrestrial atmosphere is e n = 115 GeV, while = 850 GeV. In astro- 
physical settings, the density is usually low enough so that decay always 
dominates over hadronic interactions. An intermediate case of some in- 
terest is production of secondary cosmic rays in the solar chromosphere, 
where the scale height is larger than in the Earth's atmosphere so that 
decay remains dominant for another order of magnitude (Seckel, Stanev 
& Gaisser, 1991). 

The same set of cascade equations (see Eq. 19) governs air showers 
and uncorrelated fluxes of particles in the atmosphere. The boundary 
condition for an air shower initiated by a primary of mass A and total 
energy Eq is 

N(X)\ X=0 = A5{E - Eo/A) (21) 

and iV(0) = for all other particles. This approximation, in which a 
nucleus is treated as consisting of independently interacting nucleons, 
is called the superposition approximation. In practice in Monte Carlo 
solutions of the cascade equation it is straightforward to remove this 
approximation given a model of nuclear fragmentation, (e.g. Battistoni, 
et al, 1997). 
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Figure 2. The flux of nucleons. The heavy black line shows the numerical form 
of Eq. 22. The lighter lines show extrapolations of fits (Gaisser, Honda, Lipari & 
Stanev, 2001) to measurements of protons, helium and three heavier groups below 
100 GeV/nucleon. 



The other important boundary condition is that for uncorrelated fluxes 
in the atmosphere: 

N(X)\ X=Q = <f> p (E) = 1.7 x 10 4 E- 2 - 7 (GeVcm 2 ssr) -1 . (22) 

The numerical approximation is for the flux of all nucleons summed 
over the five major nuclear groups shown in Fig. 2. In Eq. 22, E is 
total energy per nucleon. This numerical approximation is shown as 
the heavy solid line in Fig. 2. Its curvature at low energy is just a 
consequence of plotting the power law in total energy per nucleon as a 
function of kinetic energy per nucleon. Only a subset of available data is 
shown in Fig. 2. Data from the magnetic spectrometers BESS98 (Sanuki, 
et al, 2000) and AMS (Alcarez et al, 2000a) are indistinguishable 
on the plot for protons, although they differ somewhat for helium ( 
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Alcarez et al., 2000b). Data from the CAPRICE spectrometer (Boezio 
et al, 1999) are 15-20% lower than BESS98 above 10 GeV/nucleon. The 
higher energy data are from balloon-borne emulsion chambers, which are 
subject to larger systematic errors because not all energy is sampled in 
the calorimeter. Proton and helium data from JACEE (Asakimori et 
al., 1998) are shown. Data of RUNJOB (Apanasenko et al, 2001; 
Furukawa et al, 2003) are shown for five groups of nuclei (protons and 
helium, CNO, Ne-Si and Fe) above 1000 GeV. The fits to CNO, Ne-Si 
and Fe are normalized at 10.6 GeV/nucleon to data of Engelmann et 
al, 1990. 
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Figure 3. The flux of protons. The data are from the measurements of Mocchiutti, 
2003 with the CAPRICE detector. 



2.2.1 Uncorrelated fluxes in the atmosphere. The simplest 
physical example to illustrate the solution of Eq. 22 is to calculate the 
vertical spectrum of nucleons as a function of depth in the atmosphere. 



Outstanding Problems in Particle Astrophysics* 



13 



Nucleons are stable compared to the transit time through the atmo- 
sphere, so only losses due to interactions are important in the cascade 
equation 19. In the approximation of scaling, the dimensionless dis- 
tribution F N ^ N (E,E') -> F(£) as in Eq. 9, with £ = E/E' . Eq. 19 
becomes 

dN(E,X) _ N(E,X) 1 f 1 AT/ E ^ p^d^ 



- I' N(l,X)F(0§. (23) 



dA A A 

The dependence on energy and depth can be factorized, and the solution 
of Eq. 23 is 

N(E,X) = Ke~ x/AN x E~ a , (24) 

where the boundary condition 22 is satisfied if K = 1.7x 10 4 and a = 2.7. 
The attenuation length is related to the interaction length A by 

An = -r^^, (25) 

where Z^n = £ a ~ 2 Fn-*n(€) d£ ~ 0.3 is the spectrum-weighted 
moment for nucleons, analogous to Eq. 12 for pions. 

The solution outlined above has several obvious approximations such 
as neglect of nucleon-anti-nucleon production and neglect of energy- 
dependence of the cross sections, but it nevertheless gives a reasonable 
representation of measurements of the spectrum of protons at various 
atmospheric depths, as shown in Fig. 3. For comparison with measure- 
ments of protons, the solution of Eq. 24 for all nucleons must be modified 
to remove neutrons, which increase slightly as a fraction of the total flux 
of nucleons with increasing depth in the atmosphere. The correction, 
as described in Gaisser, 1990, is included in the calculations shown in 
Fig. 3. 

Another benchmark measurement of secondary cosmic rays in the 
atmosphere is the flux of muons. The main source of muons is from 
decay of charged pions. There is also a small contribution from decay 
of charged kaons, which becomes somewhat more important at high 
energy. At very high energy the muon energy spectrum becomes one 
power steeper than the parent spectrum of nucleons as a consequence of 
the extra power of E n in the ratio X^/d^, which represents the decreasing 
probability of decay relative to re-interaction for charged pions at high 
energy (see Eqs. 19 and 20). For E < e w essentially all pions decay, 
and the muon production spectrum has the same power behavior as 
the parent pion and grandparent nucleon spectrum (a ~ 2.7). At low 
energy, however, muon energy-loss and decay become important, and 
the muon spectrum at the ground falls increasingly below the production 
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Figure 4- Summary of measurements of the vertical muon intensity at the ground. 
The solid line shows an analytic calculation Gaisser, 2002. The dotted line shows the 
spectrum in the absence of decay and energy loss, or equivalently the muon production 
spectrum integrated over the atmosphere. 

spectrum. To account for all the complications one generally resorts to 
Monte Carlo calculations. However, analytic approximations (Lipari, 
1993) of these effects are also possible. The full line in fig. 4 shows one 
such calculation (Gaisser, 2002), which uses as input the simple power- 
law primary spectrum of Eq. 22. This simple result compares relatively 
well to various measurements in several energy ranges. The recent data 
are from CAPRICE (Kremer et al, 1999), L3-Cosmic (Achard et al., 
2004) and BESS (Motoki et al, 2003). One can also see the level of 
relative systematic uncertainties between different measurements. 

Although I will not discuss the subject here, The most important 
secondary cosmic-ray flux is the atmospheric neutrino beam because of 
the discovery of neutrino oscillations by Super-Kamiokande (Fukuda et 
al., 1998). The experimental situation is reviewed by Kajita & Totsuka, 
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2001 and Jung et al., 2001, and the calculations by Gaisser & Honda, 
2002. 

2.2.2 Air Showers. Above about 100 TeV cascades generated 
by individual primary nuclei have a big enough footprint deep in the at- 
mosphere to trigger an array of widely spaced detectors on the ground. 
The threshold energy may be somewhat lower for closely spaced detec- 
tors, especially at high altitude. The threshold also may be made much 
higher by separating the detectors by a large distance. Examples of the 
latter are the Akeno Giant Air Shower Array (AG AS A) (Takeda et al., 

2003) and the surface detector of the Auger Project (Watson et al., 

2004) . Such ground arrays work by looking for coincidences in an ap- 
propriate time window, then reconstructing the primary direction and 
energy from the timing pattern of the hits and the size of the signals in 
the detectors. There are large fluctuations from shower to shower which 
complicate the interpretation of the data. The air shower technique is 
used at very high energy where the flux is too low to accumulate mean- 
ingful statistics with detectors carried aloft by balloons or spacecraft. 
The dividing line at present is approximately 100 TeV. 

Because of the complicated cascade of interactions that intervenes be- 
tween the primary cosmic-ray nucleus incident on the atmosphere and 
the sparse data on the ground, Monte Carlo simulations are used to 
interpret the data. The other important reason for the necessity of a 
Monte Carlo generation of showers is that the detectors only sample a 
tiny fraction of the particles in the shower. Simulation of the response 
of an air shower detector to showers is therefore crucial. The stan- 
dard, fully stochastic, four-dimensional air shower generator is COR- 
SIKA ( Heck & Knapp, 2003). A cascade generator of similar scope and 
design is AIRES (Sciutto, 2001). A fast, one-dimensional cascade gen- 
erator (Alvarez-Muhiz et al, 2002) that uses libraries of pre-generated 
subshowers at intermediate energies inside cascades is useful for anal- 
ysis of ultra-high energy showers measured by fluorescence detectors 
for which knowledge of the lateral distribution is less important. The 
three-dimensional hybrid generator SENECA (Drescher & Farrar, 2003) 
uses stochastic Monte Carlo methods for the high-energy part of the 
shower and at the detector level, but saves time by numerically inte- 
grating the cascade equations for intermediate energies. The FLUKA 
program ( Fasso et al, 2001) is a general code for transport and inter- 
action of particles through detectors of various types, including a layered 
representation of the atmosphere. The FLUKA interaction model (Bat- 
tistoni et al, 2004) is built into the code and cannot be replaced by a 
different event generator. 
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There is a variety of hadronic event generators on the market (Bopp 
et al, 2004, Kalmykov et al, 1997, Engel et al, 2001, Werner, 1993, 
Bossard et al, 2001), which can be called by cascade programs like 
CORSIKA or AIRES to generate showers. Because the event generators 
are based on interpolations between measurements with accelerators at 
specific points in phase space and because the energies involved require 
extrapolations several orders of magnitude beyond those accessible with 
present accelerators, different hadronic event generators give different 
results for observables in air showers. We will see examples of this in 
the discussion of air showers below. 

An air shower detector essentially uses the atmosphere as a calorime- 
ter. Each shower dissipates a large fraction of its energy as it passes 
through the atmosphere, which is sampled in some way by the detector. 
It is therefore customary to plot the energy spectrum in the air shower 
regime by total energy per particle rather than by energy per nucleon 
as at low energy. In the lower energy regime the identity of each pri- 
mary nucleus can be determined as it passes through the detector on a 
balloon or spacecraft. With air showers one has to depend on Monte 
Carlo simulations to relate what is measured to the primary energy and 
to the mass of the primary particle. The resulting energy assignments 
typically have uncertainties AE / E ~ 20-30%. Primary mass is often 
quoted as an average value for a sample of events in each energy bin, or 
at best as a relative fraction of a small number of groups of elements. 

3. Acceleration 

A detailed review of the theory of particle acceleration by astrophysi- 
cal shocks is given in the lectures of Ostrowski, 2005. The main feature 
necessary for understanding the implications of air shower data for origin 
of high-energy cosmic rays is the concept of maximum energy. Diffusive, 
first-order shock acceleration works by virtue of the fact that particles 
gain an amount of energy AE oc E at each cycle, where a cycle con- 
sists of a particle passing from the upstream (unshocked) region to the 
downstream region and back. At each cycle, there is a probability that 
the particle is lost downstream and does not return to the shock. Higher 
energy particles are those that remain longer in the vicinity of the shock 
and have time to achieve high energy. 

After a time T the maximum energy achieved is 

E max ~ Zefc xBxTV s , (26) 

where f3 s = V s /c refers to the velocity of the shock. This result is 
an upper limit in that it assumes a minimal diffusion length equal to 
the gyroradius of a particle of charge Ze in the magnetic fields behind 
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and ahead of the shock. Using numbers typical of Type II supernovae 
exploding in the average interstellar medium gives E max 100 TeV 

( Lagage & Cesarsky, 1983). More recent estimates give a maximum 
energy larger by as much as an order of magnitude or more for some 
types of supernovae (Berezhko, 1996). 

The nuclear charge, Z, appears in Eq. 26 because acceleration depends 
on the interaction of the particles being accelerated with the moving 
magnetic fields. Particles with the same gyroradius behave in the same 
way. Thus the appropriate variable to characterize acceleration is mag- 
netic rigidity, R = pc/Ze ~ Et t/Ze, where p is the total momentum 
of the particle. Diffusive propagation also depends on magnetic fields 
and hence on rigidity. For both acceleration and propagation, therefore, 
if there is a feature characterized by a critical rigidity, R* , then the 
corresponding critical energy per particle is E* = Z x R*. 




E (eV/nucleus) 



Figure 5. High-energy cosmic-ray spectrum. References to the data are given in ( 
Gaisser & Stanev, 2004). The shaded region indicates a factor of 30 in total energy 
(see text). 
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4. The anatomy of the cosmic-ray spectrum 

The knee of the spectrum is the steepening that occurs above 10 15 eV, 
as shown in Fig. 5, while the ankle is the hardening around 3 x 10 18 eV. 
One possibility is that the knee is associated with the upper limit of 
acceleration by galactic supernovae, while the ankle is associated with 
the onset of an extragalactic population that is less intense but has a 
harder spectrum that dominates at sufficiently high energy. The general 
idea that the knee may signal the end of the population of particles 
produced in the Galaxy is an old one that I will trace in the next two 
subsections. 

4.1 The knee 

If the knee is a consequence of galactic cosmic accelerators reaching 
their limiting energy, then there are consequences for energy-dependence 
of the composition that can be used to check the idea. This follows 
from the form of Eq. 26. Consider first the simplest case in which all 
galactic accelerators are identical. Then E max = Z x R*, where R* 
characterizes the maximum rigidity. When the particles are classified 
by total energy per nucleus, protons will cut off first at E max = eR*, 
helium at E max = 2eR*, etc. Peters, 1961 described this cycle of 
composition change and pointed out the consequences for composition 
in a plot like that reproduced here as Fig. 6. Since the abundant elements 
from protons to the iron group cover a factor of 30 in Z, the "Peters 
cycle" should occupy a similar range of total energy. 

Because the observed spectrum does not abruptly stop, Peters hy- 
pothesized a new population of particles coming in with a somewhat 
harder spectrum, as indicated by the line in Fig. 6. Comparison with 
the data in Fig. 5 shows that reality is more complicated. The shaded 
area indicates the factor of 30 for a Peters cycle assuming R* = 10 15 eV. 
What is observed is that the steepened spectrum above the knee contin- 
ues on smoothly for more than an order of magnitude in energy before 
any sign of a hardening. Even postulating a significant contribution from 
elements heavier than iron (up to uranium, Horandel, 2004) cannot ex- 
plain the smooth continuation all the way up to the ankle. 

One possibility is that most galactic accelerators cut off around a rigid- 
ity of perhaps 10 15 eV, but a few accelerate particles to much higher en- 
ergy and account for the region between the knee and the ankle (Erlykin 
& Wolfendale, 2001). This scenario would be a generalization of Peters' 
model. Its signature would be a sequence of composition cycles alter- 
nating between light and heavy dominance as the different components 
from each source cut off. As emphasized by Axford, 1994, however, 
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Figure 6. Recreation of a figure from Peters, 1961 showing a sequence of cutoffs for 
ions with successively higher charge when classified by total energy per particle. The 
line shows a hypothetical second component with a harder spectrum than the lower 
energy component. (See text for discussion). 

the problem with this type of model is that it requires a fine-tuning 
of the high-energy spectra so that they rise to join smoothly at the 
knee then steepen to fit the data to ~ 10 18 eV. As a consequence, sev- 
eral models have been proposed in which the lower-energy accelerators 
(E < 10 15 eV) inject seed particles into another process that acceler- 
ates them to higher energy. In this way the spectrum above the knee is 
naturally continuous with the lower energy region. One such possibility 
is acceleration by interaction with multiple supernova shocks in a cluster 
of supernovae (Axford, 1994). Another possibility (Volk & Zirakashvili, 
2003) is acceleration by a termination shock in the galactic wind ( Jokipii 
& Morrill, 1987). 

The fine tuning problem (i.e. to achieve a smooth spectrum with a 
sequence of sources with different maxima) was actually clearly recog- 
nized by Peters, 1961 in his original statement of this idea. He correctly 
pointed out, however, that since the cutoff is a function of rigidity while 
the events are classified by a quantity close (but not equal) to total 
energy, the underlying discontinuities are smoothed out to some extent. 
An interesting question to ask in this context is what power at the source 
would be required to fill in the spectrum from the knee to the ankle. The 
answer depends on what is assumed for the spectrum of the sources and 
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the energy dependence of propagation in this energy region. Reasonable 
assumptions (e.g. Q(E) oc E~ 2 and r esc oc E~ s with 5 ~ 0.3) lead to an 
estimate of ~ 2 x 10 39 erg/s, less than 10% of the total power requirement 
for all galactic cosmic-rays. For comparison, the micro-quasar SS433 at 
3 kpc distance has a jet power estimated as 10 39 erg/s (Distefano et ai, 

2002) . 

Another possibility is that the steepening of the spectrum at the knee 
is a result of a change in properties of diffusion in the interstellar medium 
such that above a certain critical rigidity the characteristic propagation 
time t esc decreases more rapidly with energy. If the underlying ac- 
celeration process were featureless, then the relative composition as a 
function of total energy per particle would change smoothly, with the 
proton spectrum steepening first by 0.3, followed by successively heavier 
nuclei. It is interesting that this possibility was also explicitly recognized 
by Peters, 1961. 

A good understanding of the composition would go a long way to- 
ward clarifying what is going on in the knee region and beyond. A 
recent summary of direct measurements of various nuclei shows no sign 
of a rigidity-dependent composition change up to the highest energies ac- 
cessible (~ 10 14 eV/nucleus) (Battiston, 2004). The change associated 
with the knee is in the air shower regime. Because of the indirect nature 
of EAS measurements, however, the composition is difficult to determine 
unambiguously. The composition has to be determined from measure- 
ments of ratios of different components of air showers at the ground. 
For example, a heavy nucleus like iron generates a shower with a higher 
ratio of muons to electrons than a proton shower of the same energy. 
Swordy et al., 2002 reviewed all available measurements of the compo- 
sition at the knee. A plot of mean log mass ((\n(A))) (Ahrens et al, 
2004b) shows no clear pattern when all results are plotted together. The 
best indication at present comes from the Kascade experiment ( Roth, 

2003) , which shows clear evidence for a "Peters cycle", the systematic 
steepening first of hydrogen, then of helium, then CNO and finally the 
iron group. The transition occurs over an energy range from approxi- 
mately 10 15 eV to 3 x 10 16 eV, as expected, but the experiment runs 
out of statistics by 10 17 eV, so the data do not yet discriminate among 
the various possibilities for explaining the spectrum between the knee 
and the ankle. 

4.2 The ankle 

Above some sufficiently high energy it seems likely that the cosmic 
rays will be of extra-galactic origin. A proton of energy 10 18 eV has a 
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gyroradius of a kiloparsec in a typical galactic magnetic field, which is 
larger than the thickness of the disk of the Galaxy. Given constraints 
from observed isotropy of particles with E ~ 10 19 eV, where the cor- 
responding proton gyroradius is comparable to the full extent of the 
Galaxy, the usual assumption is that particles above 10 19 originate out- 
side our galaxy. There is a suggestion of an anisotropy just around 
10 18 eV from the central regions of the galaxy (Hayashida et al, 1999) 
that may be due to neutrons, which survive for a mean pathlength of 
~ 10 kpc at this energy, and could therefore reach us from the galactic 
center. If so, this would suggest that at least some fraction of the cosmic 
rays around 10 18 eV are still of galactic origin. An interesting discussion 
of possible implications is given by Anchordoqui et al, 2004a. In any 
case, the questions of whether there is a transition from galactic to ex- 
tragalactic cosmic rays and at what energy such a transition occurs are 
of great interest. 
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Figure 7. Four measurements of the spectrum above 10 eV . 



As background for a discussion of the significance of the ankle and 
how it should be interpreted, it is helpful to look at the spectrum from 
individual measurements separately so that systematic differences be- 
tween measurements do not obscure any details that may be present. 
Results from four groups are shown in Fig. 7: Fly's Eye stereo (Bird, 
et al, 1993); Akeno (Nagano et al, 1992) and AG AS A (Takeda et al, 
2003); Yakutsk (Glushkov et al, 2003); and the measurements made 
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with Hi-Res (Abassi et al., 2004 and Abu-Zayyad et al., 2001). In 
all these measurements there is a suggestion of a steepening just below 

10 18 eV, which is sometimes referred to as the "second knee." The ankle 
appears as a saddle-like shape with its low point between 3 x 10 18 and 

10 19 eV, depending on the experiment. One could fit the saddle with a 
final Peters cycle starting just below 10 18 eV and an extragalactic com- 
ponent crossing as in Fig. 6 to contribute to the ankle in the overall 
spectrum (Bahcall & Waxman, 2003). Alternatively, it is possible to 
make a model in which extragalactic cosmic rays account for the entire 
observed flux down to 10 18 eV (Berezinsky et al., 2004) or even lower 
( Bergman, 2004). The difference lies in the assumptions made for the 
spectrum and cosmological evolution of the sources. I will return to this 
issue in the next section. 

First it is interesting to ask what data on primary composition may 
tell us about the changing populations of particles above 10 17 eV. In this 
energy range the composition is measured by the energy dependence of 
the position of shower maximum, X max . An air shower consists of a su- 
perposition of electromagnetic cascades initiated by photons from decay 
of 7T particles produced by hadronic interactions along the core of the 
shower as it passes through the atmosphere. Most of the energy of the 
shower is dissipated by ionization losses of the low-energy electrons and 
positrons in these subshowers. The composite shower reaches a max- 
imum number of particles (typically 0.7 particles per GeV of primary 
energy) and then decreases as the individual photons fall below the crit- 
ical energy for pair production. Because each nucleus of mass A and 
total energy Eq essentially generates A subshowers each of energy Eq/A 
the depth of maximum depends on Eq/A. Since cascade penetration 
increases logarithmically with energy, 

X max = \ ER \og(E Q /A) + C, (27) 

where Xer is a parameter (the "elongation rate") that depends on the 
underlying properties of hadronic interactions in the cascade. 

Fig. 8 shows results of measurements with the Fly's Eye stereo de- 
tector (Bird, et al., 1993) compared to measurements of HiRes (Abu- 
Zayyad et al., 2001; Archbold & Sokolsky, 2003). A weak inference 
about composition can be made by comparing to the results of simu- 
lations, two of which are shown in the figure. Both calculations use 
CORSIKA ( Heck k, Knapp, 2003) with two different interaction mod- 
els (Kalmykov et al., 1997, Engel et al., 2001). The measurement with 
the Fly's Eye Stereo detector suggests a transition from a large fraction 
of heavies below 10 18 eV to a larger fraction of protons by 10 19 eV (how 
much larger depending on which interaction model is chosen) . The coin- 
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Figure 8. Plot of data on mean depth of maximum vs energy. Filled squares 
are data from the stereo Fly's Eye (Bird, et ai, 1993). Open symbols show the 
data of HiRes (squares) (Archbold & Sokolsky, 2003) and HiRes prototype with MIA 
(circles) (Abu-Zayyad et ai, 2001). This figure also appears in Ref. (Gaisser & 
Stanev, 2005). 



cidence of the change of composition from heavier toward lighter) around 
~ 3 x 10 18 eV with the ankle feature in the Fly's Eye data at the same 
energy led to the suggestion of a transition from galactic to extragalactic 
cosmic rays as a possible explanation (Bird, et al., 1993). This interpre- 
tation would favor a model like that of Bahcall & Waxman, 2003. The 
more recent HiRes data set, however, shows the transition from heavier 
toward lighter beginning at 10 17 eV and complete by 10 18 eV, consistent 
with the models of Berezinsky et al, 2004 and Bergman, 2004. 

Because of the uncertainties in the interaction models above acceler- 
ator energies coupled with statistical and systematic limitations of the 
experiments, the primary composition as a function of energy above 
10 17 eV remains an open question (Watson, 2004). 
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5. Highest energy cosmic rays 

Protons lose energy by three loss processes during propagation through 
the cosmos. Red-shift losses (adiabatic losses due to expansion of the 
Universe), which apply to all particles, become important when the dis- 
tance scales are comparable to the Hubble distance f» 4Gpc. Protons of 
sufficiently high energy also interact with the microwave background and 
lose energy to electron-positron pair production (for E > ~ 10 18 eV) 
and to photopion production (for E >~ 5 x 10 19 eV). The corre- 
sponding attenuation lengths (for reducing energy by a factor 1/e) are 
A e + e - ~ 1 Gpc and ~ 15 Mpc, respectively. The photo-pion process 
leads to the expectation of a suppression of the flux above 5 x 10 19 eV 
unless the sources are within a few tens of Mpc. The suppression is 
referred to as the GZK cutoff (or GZK feature) in recognition of the 
authors of the two papers, Greisen, 1966 and Zatsepin & Kuz'min, 1966 
who first pointed out the effect shortly after the discovery of the mi- 
crowave background radiation. 

The actual shape of the spectrum at Earth after accounting for these 
three loss processes depends on what is assumed 

■ for the spatial distribution of sources, 

■ for the spectrum of accelerated particles at the sources, and 

■ for the possible evolution of activity of the sources on cosmological 
time scales. 

A classic calculation is that of Berezinsky Sz Grigorieva, 1988, in which 
the energy-loss equation is integrated numerically. This approach ne- 
glects effects of fluctuations, which may be noticeable in certain circum- 
stances. A recent example of a Monte Carlo propagation calculation 
of cosmological propagation is by Stanev et ai, 2000, which contains 
comparisons with other calculations. Figure 9 from Abassi et ai, 2002 
shows an example of a calculated cosmologically evolved spectrum com- 
pared to data of HiRes (Abassi et ai, 2004) and AGASA (Takeda et 
ai, 2003). 

The calculation illustrated in Fig. 9 is for a uniform distribution of 
sources out to large red shifts. The "pile-up" (which is amplified by 
plotting E 3 x differential flux) is populated by particles that have fallen 
below the GZK energy. The saddle or ankle feature on the E 3 dN / dE 
plot is a consequence of energy losses to pair production. The data be- 
low 10 18 eV have been fit in this example by adding an assumed galactic 
contribution, which is not shown separately here (but see Bergman, 2004 
for a discussion). As pointed out long ago by Hillas, 1968 (see also Hillas, 
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Figure 9. The ultra-high-energy cosmic-ray spectrum from a paper by Abassi et 
ai, 2002. The points with the fitted curve are from the HiRes fluorescence detectors ( 
Abassi et ai, 2004), while the higher set of points are the AGASA spectrum (Takeda 
et at, 2003). See text for a discussion of the curve. 

1974), the degree to which the extragalactic flux contributes at low en- 
ergy (for example below 10 18 eV) can be adjusted by making different 
assumptions about increased activity at large redshift, z > 1. Such high 
redshift sources do not contribute at high energy because of adiabatic 
losses, but they can contribute at low energy if they were sufficiently 
abundant. On the other hand, it is also possible to construct a model 
in which the extragalactic sources only contribute significantly at very 
high energy (e.g. above 10 19 eV), as in Bahcall & Waxman, 2003. 

The question of the energy above which the extragalactic component 
accounts for all of the observed spectrum is important because it is 
directly related to the amount of power needed to supply the extra- 
galactic component. An estimate of the power needed to supply the 
extragalactic cosmic rays is obtained by the replacements r esc — > th and 
(f)(E) — > 4>eg(E) in Eq. 6, where th ~ 1-4 x 10 10 yrs is the Hubble 
time and (Peg is the extragalactic component of the total cosmic-ray 
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flux. The result depends on where the EG-component is normalized to 
the total flux and on how it is extrapolated to low energy. A minimal 
estimate follows from taking the hardest likely spectral shape (differen- 
tial index a = —2) and normalizing at 10 19 eV. This leads to a total 
power requirement of Peg = 10 37 erg/Mpc 3 /s. If the rate of GRBs 
is 1000 per year, then they would need to produce 2 x 10 52 erg/burst 
in accelerated cosmic rays to satisfy this requirement. The density of 
AGNs is ~ 10~ 7 /Mpc 3 (Peebles, 1993), so the corresponding require- 
ment would be 10 44 erg/s per AGN in accelerated particles. Normalizing 
the extragalactic component at 10 18 eV increases the power requirement 
by a factor of 10, and assuming a softer spectral index increases it even 
more. These are very rough estimates just meant to illustrate the re- 
lation between the intensity of extragalactic cosmic rays and the power 
needed to produce them. A proper treatment requires an integration 
over redshift, accounting for the spatial distribution of sources and their 
evolution over cosmic time (Stanev, 2003). 

In addition to being able to satisfy the power requirement, the sources 
also have to be able to accelerate particles to ~ 10 20 eV. The require- 
ments on size and magnetic field in the sources are essentially given by 
setting -Emax = 10 20 eV in Eq. 26 and solving for the product of B x T V s 
(or B x R). The resulting constraint on the sources is the famous plot 
of Hillas, 1984. Sites with a sufficiently large B x R included active 
galactic nuclei (AGN) (Berezinsky et ai, 2004) and the termination 
shocks of giant radio galaxies (Biermann & Strittmatter, 1987). Since 
1984 two more potential candidates have emerged, jets of gamma-ray 
bursts (GRB Waxman, 1995, Vietri, 1995) and magnetars. 

Neutrino telescopes may make a decisive contribution to identifying 
ultra-high energy cosmic accelerators by finding point sources of TeV 
neutrinos and associating them with known objects, for example with 
gamma-ray sources such as AGN or GRB. Even the detection of a dif- 
fuse flux of high-energy, extraterrestrial neutrinos, or the determination 
of good limits can constrain models of cosmic-ray origin under the as- 
sumption that comparable amounts of energy go into cosmic-rays and 
into neutrinos produced by the cosmic-rays at the sources (Anchordoqui 
et al, 2004b). 

As is well-known, the AGASA measurement (Takeda et al., 2003) 
shows the spectrum continuing beyond the GZK energy without a sup- 
pression, while the HiRes measurement (Abassi et al, 2004) is consistent 
with the GZK shape. The discrepancy between the two experiments is 
not as great as it appears in Fig. 9. The difference is amplified quadrat- 
ically when the differential spectrum is multiplied by E 3 in the plot. A 
20-30% shift in energy assignment (downward for AGASA or upward 
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for HiRes) brings the two data sets into agreement below 5 x 10 19 eV. 
A quantitative estimate of the statistical significance of the difference 
above 5 x 10 19 eV is ambiguous because it depends on what is assumed 
for the true spectrum. The difference is generally considered to be less 
than three sigma, however (Olinto, 2004). One of the first goals of the 
Auger project (Klages, 2005) is to resolve the question of whether the 
spectrum extends significantly beyond the GZK energy. If not, then 
higher statistics accumulated by Auger may be used to constrain mod- 
els of extra-galactic sources by making precise measurements both above 
and below the GZK energy. If the spectrum does continue beyond the 
cutoff then identification of specific, cosmologically nearby sources from 
the directions of the ultra-high energy events would be likely. 

6. Conclusion and outlook 

While the questions surrounding the highest energy particles are clearly 
of the greatest importance, there are several important open questions 
at all energies. One is the transition, if it exists, from galactic to extra- 
galactic cosmic rays. Proposed experiments associated with the tele- 
scope array (Arai et al., 2003) may be optimized to provide good cov- 
erage down to 10 17 eV (Thompson, 2004). 

Understanding the knee of the spectrum remains an outstanding prob- 
lem in cosmic-ray astrophysics. Kascade-Grande (Huangs et al, 2003) 
will cover the energy range from below the knee to 10 18 eV with a multi- 
component air shower array at sea level. The IceCube detector at the 
South Pole (Ahrens et al, 2004a) will have a kilometer square surface 
component, IceTop (Gaisser et al, 2003; Stanev, et al, 2005), 1.4 
km above the top of the cubic kilometer neutrino telescope. The whole 
constitutes a novel, large three-dimensional air shower detector with 
coverage of the cosmic-ray spectrum from below the knee to 10 18 eV. 
In addition to its main mission of neutrino astronomy, IceCube there- 
fore also has the potential to make important contributions to related 
cosmic-ray physics. The high altitude of the surface (9300 m.a.s.l.) al- 
lows good primary energy resolution with the possibility of determining 
relative importance of the primary mass groups from the ratio of size 
at the surface to muon signal in the deep detector. With a coverage 
extending to 10 18 eV, these detectors also have the potential to clarify 
the location of a transition to cosmic-rays of extragalactic origin. 

While the knee of the spectrum remains in the realm of air shower 
experiments for the time being, there is much activity aimed at ex- 
tending direct measurements of the primary spectrum and composi- 
tion to reach the knee from below. The detectors include Advanced 
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Thin Ionization Calorimeter (ATIC, Wefel et al, 2003); Transition 
Radiation Array for Cosmic Energetic Radiation (TRACER, Midler, 
2005); Cosmic Ray Energetics and Mass (CREAM, Seo et al, 2003) 
and Advanced Cosmic-ray Composition Experiment for Space Station 
(ACCESS, http://www.atic.umd.edu/access.html). While ACCESS is 
to be flown in space, ATIC, TRACER and CREAM all take advantage 
of NASA's long-duration balloon program which is regularly achieving 
flights of two to four weeks in Antarctica. These experiments also have 
the opportunity to extend measurements of the ratio of secondary to 
primary nuclei to much higher energy and hence to resolve the related 
questions about the average source spectral index and about the isotropy 
of galactic cosmic rays. 
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